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Universal quantum computation can be achieved by simply performing single-qubit measurements 
on a highly entangled resource state, such as cluster states. The family of AfHeck-Kennedy-Lieb- 
Tasaki (AKLT) states has recently been intensively explored and shown to provide restricted com- 
putation. Here, we show that the two-dimensional AKLT state on a honeycomb lattice is a universal 
resource for measurement-based quantum computation. 



PACS numbers: 03.67.Lx, 75.10.Jm, 64.60.ah 

Introduction. Quantum computation promises expo- 
nential speedup over classical computation by exploiting 
the quantum mechanical nature of physical processes [T] . 
In addition to the standard circuit model based on uni- 
tary evolution, surprisingly, local measurement alone 
provides the same power of computation, given only a 
prior sufficiently entangled state [2] [3] . For this model 
of measurement-based quantum computation (MBQC), 
universal resource states are known to be very rare [4], 
but examples do exist 5-8 . The 2D cluster state on 
the square lattice is a universal resource state [21 [S]. 
Cluster states can be created by the Ising interaction 
from uncntangled states [9], but they do not arise 
as unique ground states of two-body interacting Hamil- 
tonians [TU]- However, by careful design of Hamilto- 
nians, certain ground states can be used for universal 
MBQC [lUIIl], and this opens up an appealing possibil- 
ity of creating universal resource states by cooling. 

A new perspective on MBQC emerged when it was dis- 
covered that the one-dimensional Affleck-Kennedy-Lieb- 
Tasaki (AKLT) state [13], originally constructed in the 
setting of condensed matter physics, can serve as re- 
sources for restricted computations [6] [14] [15] . In any di- 
mension, the AKLT state is the ground state of a particu- 
larly simple Hamiltonian which only has nearest-neighbor 
two-body interactions, is rotationally invariant in spin 
space and shares all spatial symmetries of the underlying 
lattice [T5] . The discovery of the resourcefulness of AKLT 
states creates additional avenues for its experimental re- 
alization [17] . and has instilled novel concepts in MBQC, 
such as the renormalization group and the holographic 
principle [T5J[TH]. However, to date one crucial element 
was missing in this direction: the AKLT family was not 
known to contain a universal resource. Here, we over- 
come this gap by demonstrating that the AKLT state 
on a two-dimensional honeycomb lattice is a universal 
resource for measurement-based quantum computation. 

To do this, we proceed in three steps. First, we show 
that it can be mapped to a random planar graph state 
\G(A)) by local generalized measurement, with the graph 
G(A) depending on the set A of measurement outcomes 
on all sites (defined below). Second, we argue that the 
computational universality of a typical resulting graph 
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FIG. 1: Illustrations of the AKLT state and the 2D cluster 
state, (a) AKLT state. Spin singlets of two virtual spins 1/2 
(i.e. qubits) are located on the edges of the honeycomb lattice. 
A projection Ps, v at each lattice site v onto the symmetric 
subspace of three virtual spins defines the AKLT state. In 
one hexagon, sites are labeled by A or B to show the bi- 
partite (or bi-colorable) property of the honeycomb lattice, 
(b) 2D Cluster state. One qubit resides at each lattice site 
and one stabilizer generator is shown. 

state |G(„4)) hinges solely on the connectivity of G(A), 
and is thus a percolation problem. Third, we demon- 
strate via Monte Carlo simulation that the typical graphs 
G(A) are indeed deep in the supercritical phase. 

The AKLT state [T3J on the honeycomb lattice C has 
one spin-3/2 per site of C. The state space of each spin 
3/2 can be viewed as the symmetric subspace of three 
virtual spin-l/2's, i.e., qubits. In terms of these virtual 
qubits, the AKLT state on C is (see Fig. [iji) 

|$AKLT>= (g) P S ,v <S> W«> C 1 ) 
veV(C) e€E(C) 

where V{C) and E(C) to denote the set of vertices and 
edges of C, respectively. Ps, v is the projection onto the 
symmetric (equivalently, spin 3/2) subspace at site v of 
C [20]. For an edge e = (v,w), \4>) e denotes a singlet 
state, with one spin 1/2 at vertex v and the other at w. 

A graph state \G) is a stabilizer state [21] with one 
qubit per vertex of the graph G and is the unique eigen- 
state of a set of commuting operators [5], usually called 
the stabilizer generators [2"2"| . 
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TABLE I: The dependence of stabilizers and encodings on the 
local POVM outcome. \C\ denotes the total number of sites 
contained in a domain and i,j = 1..3jCj. The honeycomb 
lattice C is bi-partite and all sites can be divided into either 
A or B sublattice, V(C) = Au B; see Fig. One choice of 
the sign is A; — 1 if the virtual qubit i £ v £ A and Ai = — 1 
if i £ v' £ B; this is due to the negative sign in the stabilizer 
generator for a singlet, i.e., (— ct m ta^ j)\4>}ij = \4>)ij f° r an 



where nb(u) denotes the neighbors of vertex v, and 
X = o x , Y = a y and Z = o~ z are the Pauli matrices. A 
cluster state is a special case of graph states, with the un- 
derlying graph being a regular lattice (see Fig. [TJd) . Any 
2D cluster state is a universal resource for measurement- 
based quantum computation [2l I23|. 

Reduction to a graph state. To show that the 2D AKLT 
state of four-level spin-3/2 particles can be converted to 
a graph state of two-level qubits, we need to preserve 
a local two-dimensional structure at each site. This is 
achieved by a local generalized measurement [T], also 
called positive-operator- value measure (POVM), on ev- 
ery site v on C. The POVM consists of three rank-two 
elements 
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(|000)(000| + | 111) <1H|) (3a) 

(I +++>(+ + +| + I )( 1) (3b) 

(\i,i,i)(i,i,i\ + | -i, -4, -£, -^|), (3c) 



where |0/1), |±) = (|0) ± |1»/V2 and | ± i) = (|0) ± 
i|l))/\/2 are eigenstates of Z, X and Y, respectively, and 
|0) = | t}> |1) = I D- Physically, F V)0 is proportional to a 
projector onto the two-dimensional subspace spanned by 
the S a = ±3/2 states. The above POVM elements obey 
the relation J2ve{x,y,z} F v,u F v,v = Ps,v, i-e-, project onto 
the symmetric subspace, as required. The outcome a v of 
the POVM at site v is random, which can be x, y or z, and 
it is correlated with the outcomes at other sites due to 
the entanglement in the AKLT state [13l [24] . As demon- 
strated below, the resulting quantum state, dependent 
on the random POVM outcomes A — {a v , v € V(C)}, 

l*M)> = Fv,a v |*aklt) = (g) F v , av (g) \4>) e 
vev(c) vev(c) eeE(C) 

(4) 

is equivalent under local unitary transformations to an 
encoded graph state \G(A)). We show that the corre- 
sponding graph G(A) is constructed from graph C by 
applying the following two rules: 
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FIG. 2: Graphical rules for transformation of the lattice C 
into the graph G(A), depending on the POVM outcomes A. 
a) Rules Rl and R2. b) An example of obtaining G(A). The 
alphabets inside the circle indicate the POVM outcomes, c) 
An example to illustrate the qubit encoding of a domain, d) 
An example for demonstrating the stabilizer generator. 



Rl (Edge contraction): Contract all edges e G E(£) 
that connect sites with the same POVM outcome. 

R2 (Mod-2 edge deletion): In the resultant multi- 
graph, delete all edges of even multiplicity and con- 
vert all edges of odd multiplicity into conventional 
edges of multiplicity 1 . 

A set of sites in C that are contracted into a single ver- 
tex of G(A) by rule Rl is called a domain. Each domain 
supports a single encoded qubit. The stabilizer gener- 
ators and the encoded operators for the resulting codes 
are summarized in Table [I] 

Rule Rl derives intuitively from the antiferromag- 
netic property of the AKLT state: neighboring spin- 
3/2 particles must not have the same S a = 3/2 (or - 
3/2) configuration [13]. Hence, after the projection onto 
S a = ±3/2 subspace by the POVM, the configurations 
for all sites inside a domain can only be 1 3/2, —3/2, ... ) 
or | — 3/2, 3/2, . . . ), and these form the basis of a single 
qubit. This can also be understood in terms of the stabi- 
lizer. Consider the case where two neighboring POVMs 
yield the same outcome, say z; see Fig. [2];. Due to 
the projections F u z and F ViZ (with u — {1,2,3} and 
v = {4,5,6} each containing three virtual qubits), the 
operators Z\Z%, Z2Z3, and Z§Z§ become stabilizer 

generators of |^(^4)). Moreover, the stabilizer — Z3Z4 of 
|0)34 commutes with F u>z ® F VjZ , and thus remains a sta- 
bilizer element for l'F(.A)) [20]. In brief, the stabilizer 
generators Z1Z2, Z2Z3, —Z3Z4, Z4Z5, Z§Z§ lead to a sin- 
gle encoded qubit 
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supported by the two sites u and v jointly. Note the 
antiferromagnetic ordering [13 among groups of three 
virtual qubits. To reduce the support of this logical 
qubit to an individual site, a measurement in the basis 
{|(000)„) ± |(lll)t,}} is performed. The resulting state is 
a|(000)„)±/3|(lll) M ), with the sign "±" known from the 
measurement outcome. This is the proper encoding for a 
domain consisting of a single site. Domains of more than 
two sites are thereby reduced to single sites. 

To see that | , 3>(„4)) is indeed equivalent under local uni- 
tary transformations to an encoded graph state \G(A)), 
we consider the example of four domains c, u, v, w, each 
consisting of a single site of £, where the POVM out- 
come is z on the central domain c and x on all lateral 
domains u, v and w; see Fig.[2ji. By similar arguments as 
above [2D], the operator O = — X\X V X2X1' X-^X^ is in 
the stabilizer of ^(.A)). Using the encoding in Table [Ij 
i.e., with the encoded Pauli operators X c = Z\ZiZz, 
Z u = ±Xy , Z v = ±X2' , and Z w = iX^/, we find that 
O = ±X C Z U Z V Z W which is (up to a possible sign) one 
stabilizer generator defining the graph state; see Eq. ([2]). 

By the above construction, if two domains it, v are 
connected by an edge of multiplicity m, the inferred 
graph state stabilizer generators will contain factors of 
X U Z V or X V Z U . Rule R2 thus follows the observa- 
tion that Z 2 = I. Generalizing these ideas, one can rig- 
orously prove that for any A of POVM outcomes, the 
state |^(^4)) is local-unitarily equivalent to an encoded 

We shall denote by \G(A)) the 



graph state \G(A)) [25] 
same graph state but with domains of single sites. 

Random graph states and percolation. Whether or not 
typical graph states \G(A)) are universal resources hinges 
solely on connectivity properties of G(A), and is thus 
a percolation problem [26) . Specifically, for a large ini- 
tial C the random graph state \G(A)) can, with close 
to unit probability, be efficiently reduced to a large two- 
dimensional cluster state if the following properties hold: 

CI The distribution of the number of sites in a domain 
(i.e. domain size) is microscopic, i.e., the largest 
domain size can at most scale logarithmically with 
the total number of sites |V(£)| in the large C limit . 

C2 The probability of the existence of a path through 
G(A) from the left to the right (or top to bottom) 
approaches unity in the limit of large C 

Condition CI ensures that the graph G(A) remains 
macroscopic if the original C was, and Condition C2 en- 
sures that the system is in the supercritical phase with a 
macroscopic spanning cluster. 

Together with planarity, which holds for all graphs 
G(A) by construction, the conditions CI and C2 are suf- 
ficient for the reduction of the random graph state to a 
standard universal cluster state. The proof [35] extends 
a similar result already established for site percolation on 
a square lattice |27j . The physical intuition comes from 
percolation theory. In the supercritical phase (where 
there exists a macroscopic spanning cluster and connects 



3.6 



3.2 



s 

a 2.8 



2.4 



ave domain size 
ave distr width 
ave vertex deg 



30 



20 



10 



20000 
N 



40000 



0.01 



0.02 



0.03 



0.04 



0.05 



1/L 



FIG. 3: (color online) Average domain size, average width 
of domain size distribution, average degree of a vertex, and 
the largest domain size (inset) in the typical graphs vs. L, 
with N = L 2 being the total number of sites. For better 
discernibility, we suppress the errorbars for one set of data. 
The largest domain size scales with N as 3.337 ln(iV) — 5.566. 
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FIG. 4: (color online) Percolation study of the random graphs 
of domains: probability of a spanning cluster p c luster vs. that 
to delete a vertex Pdcictc- The threshold for destroying the 
spanning cluster is around pdcietc ~ 0.33 in deleting vertices. 



one boundary to the other), the spanning cluster con- 
tains a subgraph which is topologically equivalent to a 
coarse-grained two-dimensional lattice structure. This 
subgraph can be carved out and subsequently cleaned off 
all imperfections by local Pauli measurements, leading to 
a perfect two-dimensional lattice. 

Numerical results. We used Monte Carlo simulations 
to sample typical random graphs resulting from the 
POVM and compute their properties. The simulations 
utilize a generalized Hoshen-Kopelman algorithm [25J to 
identify domains. Due to the entanglement in the AKLT 
state, the local POVM outcomes are correlated which is 
fully taken into account in our simulations. In particular, 
to sample typical POVM outcomes A correctly, we use 
a Metropolis method to update configurations. For each 
site we attempt to flip the type (either x, y or z) to one 
of the other two equally and accept the flip with a proba- 
bility p acc = min{l,2l v 'H£'Hi'1+|£|} ! wne re \V\ and \£\ 
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denote the number of domains and inter domain edges, 
respectively, before the flip and Rule 2, and similarly \V'\ 
and \£'\ for the flipped configuration [2"5] . 

We have analyzed lattices of size up to 200 x 200 sites. 
The average degree of vertices in the typical random 
graphs is about 3.52(1), when extrapolated to an infinite 
system size. Furthermore, the typical random graphs 
retain large number of vertices \V\ = 0.495(2)iV, edges 
\E\ = 0.872(4)A, and independent cycles (or the Betti 
number) B = 0.377(2)iV, where N = L x L is the total 
number of sites in the initial honeycomb lattice. The size 
of the largest domain was never macroscopic and followed 
a logarithmic dependence on N. The average number 
of sites v £ V(C) contained in a typical domain, when 
extrapolated to the infinite system, is 2.02(1) and the 
width in the domain size distribution is extrapolated to 
1.94(1). Our simulations show that condition CI holds; 
see Fig. [3j For all POVM outcomes sampled, a macro- 
scopic cluster existed, allowing a horizontal and a vertical 
traversing path through the resulting graphs G(A). This 
shows that condition C2 holds. 

Robustness. We now quantify how deep typical graphs 
G(A) are in the supercritical phase of the percolation 
transition. A first measure is the average vertex degree. 
A heuristic argument based on a branching process sug- 



gests that a graph has a macroscopic connected compo- 
nent whenever the average vertex degree is d > 2 [29j . 
In our case, the typical graphs G(A) have an average de- 
gree of 3.52, suggesting that the system is deep in the 
supercritical phase. Furthermore, we randomly delete a 
fraction of vertices or edges from G(A). On average, it re- 
quires a deletion probability as high as pJeiete = 0.33(1) 
for vertices (see Fig. [4| and Pd clctG = 0.43(1) for edges 
(not shown) in order to the spanning cluster. These 
numbers demonstrate the robustness of the connectivity 
property. 

Concluding remarks. We investigated MBQC on the 
AKLT states and established one crucial missing ingre- 
dient in this area: the two-dimensional spin-3/2 AKLT 
state on a honeycomb lattice is indeed a universal re- 
source. The approach described in this work also applies 
to other trivalent lattices, such as the Archimedean lat- 
tices: (3, 12 2 ), (4,6,12) and (4,8 2 ), which have higher 
percolation thresholds than the honeycomb lattice. 

After the completion of our work, we learned of a sim- 
ilar result by Miyake with a different approach [3D]. 
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Appendix A: Projection onto the symmetric subspace of three qubits 

. The addition of angular momenta for three qubits (spin-1/2) gives rise to three subspaces: = \ © | © §. 

The four basis states of the S = 3/2 subspaces are |3/2, 3/2), |3/2, —3/2), |3/2, 1/2) and |3/2,-l/2) (with the 
quantization axis being assumed to be z), and they can be expressed in terms of the three-qubit basis states (with 
|0/l) = |m> = |l/2,±l/2)) 

13/2,3/2) = |000), |3/2,-3/2) = |111), (Ala) 
|3/2, 1/2) = \W) ee 4(|001) + |010) + |100)), (Alb) 

|3/2, -1/2) = \W) EE ^(|H0) + |101) + |011)). (Ale) 

As can be clearly seen, states in this subspace are symmetric under permutation of the three qubits and the corre- 
sponding projection operator is 

P s , v ee |000)(000| + \W)(W\ + \W)(W\ + |111)(1H|. (A2) 

The POVM and the post-POVM state. To show that the 2D AKLT state of four-level spin-3/2 particles can be 
converted to a graph state of two- level qubits, we need to preserve a local two-dimensional structure at each site. This 
is achieved by a local generalized measurement [I], also called positive-operator- value measure (POVM), on every site 
v on C. The POVM consists of three rank-two elements F^ a F v ^ a with a — x,y, or z, and 

|(|000)(000| + I 111) <111|) (A3a) 

F v , x = +++)<+ + +| + |-—){ 1) (A3b) 

F v ,y = \J^(\i,i,i)(i,i,i\ + \-i,-i,-i)(-^i,-i,-i\), (A3c) 

where |0/1), |±) = (|0) ± |l))/\/2 and | ±i) = (|0) ± i\l))/y/2 are eigenstates of the three Pauli operators Z, X and 
Y, respectively. Physically, F v a is proportional to a projector onto the two-dimensional subspace spanned by the 
S a = ±3/2 states, i.e., 



7 «,o = y|(|3/2,+3/2) a (3/2,+3/2| + |3/2, -3/2) Q (3/2, +3/2|), 



(A4) 



with 5*013/2, ±3/2) a = ±3/2|3/2, ±3/2) a and a indicating the quantization axis. 

The above POVM elements obey the relation ^2 u£ { x y 2 j F^ v F VjV — Ps, v — ls=3/2, i-e., project onto the symmetric 
subspace, as required. The outcome a v of the POVM at site v is random, which can be x, y or z. The resulting 
quantum state, dependent on the random POVM outcomes A = {a v ,v e V(C)}, 

l*G4» = (8) F v,- I*aklt) = (g) F v>av (g) |0) e , (A5) 

v£V(C) veV(C) e£E(C) 

where \4>) e denotes a singlet state (|01) — 1 10) ) / on the edge e. 

Example 1: qubit encoding. Consider the case where two neighboring POVMs yield the same outcome, say z; see 
Fig. 2c in the main text. Due to the projections F UtZ and F v ^ z (with u — {1,2,3} and v = {4,5,6} each containing 
three virtual qubits), the operators Z1Z2, Z2Z3, and Z4Z5, Z§Zq become stabilizer generators of |W(^4)), as, e.g., 

ZiZ 2 F u , z = F u>z (A6) 

and hence 

Z1Z2MA)) =Z X Z 2 (g) F Vtav (g) \</>) e = \H((A)). (A7) 

veV(C) e£E(C) 



6 



Moreover, the stabilizer —Z 3 Z 4 of |<^) 34 (i.e., — Z 3 Z4|0) 3 4 = 1^)34) commutes with F UyZ <g> F v>z , and thus remains a 
stabilizer element for |\I>(.4)), 

-Z 3 Z 4 \*(A))= (g) F VtUv (-Z 3 Z 4 ) (g) \d>) e = \*(A)). (A8) 

vEV(C) eEE(C) 

Therefore, regarding sites u and v, |(000) U (111)„) and |(111)„(000)„) are the only two basis states that are stabilized 
(i.e., common eigenstates with eigenvalue equal to unity) by the above five operators, Z\Z 2 , Z 2 Z 3 , —Z3Z4, Z 4 Z 5 , Z 5 Zq. 
This leads to a single encoded qubit 

a|(000) u (lll)„) + /3|(111)„(000)„>, (A9) 

supported by the two sites u and v jointly. If more neighboring sites share the same POVM outcome, the encoded 
qubit can be easily extended. 

Example 2: graph-state stabilizer generator . To see that ^(A)) is indeed equivalent under local unitary transfor- 
mations to an encoded graph state \G(A)), we consider the example of four domains c, u, v, w, each consisting of a 
single site of C, where the POVM outcome is z on the central domain c and x on all lateral domains u, v and w; see 
Fig. 2d of the main text. By direct computation we have that 

[— XiXy, F UyX ] = [— X 2 X 2 >, F ViX ] = [— X 3 X 3 i , F w>x ] = 0. (A10) 

Although — XiXii individually does not commute with F CyZ , the operator O = X\XyX 2 X 2 >X 3 Xy does, as by direction 
computation, 

XiX 2 X 3 (|000)(000| + I 111) (111|) = |111)(000| + |000)(111| = (|000)(000| + |lll)(lll|)XiX 2 X 3 , (All) 
we have 



[X 1 X 2 X 3 ,F c J=0. (A12) 

This shows that O is in the stabilizer of \^(A)),_i.e., 0\$>(A)) = \^(A)}. JJsing the encoding in Table 1 of the 
main text, i.e., with the encoded Pauli operators X c = ZiZ 2 Z 3 , Z u = Xy , Z v = X 2 i, and Z w = Xy, we find that 
O = X C Z U Z V Z W which is (up to a possible sign) one stabilizer generator defining the graph state. 

The Hamiltonian and boundary conditions. The construction of the AKLT state by the valence-bond picture gives 
rise to a state that is the ground state of the following Hamiltonian 

h= e ^r 3) = e [^■^ + ^-^ 2+ ^- s ^ 3+ ml (Ai3) 

edge edge 

where is a projector of the neighboring sites i and j onto a total S = 3 subspacc. In the case of the periodic 

boundary condition, the AKLT is the unique ground state. In the case of the open boundary condition, one can 
terminate every boundary spin-3/2 by a spin-1/2, and add a corresponding term in the Hamiltonain, 

+ % (A14) 

where Si is the spin-3/2 operator at the boundary and sV is the associated spin-1/2 operator. These additional 
Hamiltonian terms for all boundary pairs make the AKLT state a unique ground state. 

In the case of the open boundary condition, the planarity of the graph is preserved even after the POVMs. However, 
in the case of the periodic boundary condition, the underlying topology is that of a torus. To make the graph of 
the corresponding graph state after the POVMs be planar, one simply measures the logical Z on sites along the two 
independent cycles and this will cut the torus into a plane. 



